Orbital order out of spin disorder: How to measure the orbital gap 
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The interplay between spin and orbital degrees of freedom in the Mott-Hubbard insulator is studied 
by considering an orbitally degenerate superexchange model. We argue that orbital order and the 
orbital excitation gap in this model are generated through the order-from-disorder mechanism known 
previously from frustrated spin models. We propose that the orbital gap should show up indirectly 
in the dynamical spin structure factor; it can therefore be measured using the conventional inelastic 
neutron scattering method. 



PACS numbers; 75.10.-b, 75.30.Ds, 75.30.Et 

The recent renaissance in the study of transition metal 
oxides has emphasized the important role being played 
by the orbital (pseudo) degeneracy inherent to perovskite 
lattices. First of all, the type of spin structure and the 
character of spin excitations [CKUcially depend on the ori- 
entation of occupied orbitalsJlTEl. Second, the excitations 
in the orbital sector get coupled to the other degrees 
of freedom (electronic, lattice, spin) and might therefore 
strongly modifKiheir excitation spectra. It was also sug- 
gested recentlyaa that low-energy orbital fluctuations are 
responsible for the highly correlated metallic state of fer- 
romagnetic manganites. 

Apparently, orbital order and orbital fluctuations de- 
serve for careful theoretical and experimental study. 
However, the orbital excitation itself is spin- and charge- 
less and can therefore be detected only indirectly due to 
its coupling to. other types of excitations. For instance, 
Ishihara et alu have recently discussed the possibility to 
detect orbital excitations by means of the anomalous x- 
ray scattering method. In this paper we propose the idea 
of detecting the orbital excitation in a conventional in- 
elastic neutron scattering experiment. That is, due to 
the inherent coupling of orbitals to the spin, a magnon 
and a single orbital wave can be excited by neutrons. Be- 
low, we illustrate this idea by considering an eg orbitally 
degenerate antiferromagnetic model. 

Physically realistic spin-orbital modelgj'El are usually 
rather complicated. To make the discussion more trans- 
parent we considjCr a simplified version of the Kugel- 
Khomskii modelffl 



1 



(1) 



describing the superexchange process in a Mott-Hubbard 
insulator with eg degeneracy, where the Hund's splitting 
in the intermediate state has been neglected. Here, Si 
is the spin 1/2 operator, while operators rf act in the 



orbital subspace with basis vectors (1,0), (0,1) corre- 
sponding to the eg{x'^ — V^) x) and eg{3z'^ — r^) ~| z) 
orbital states, respectively. The structure of r" depends 
on the index a which specifies the orientation of the bond 
(ij) relative to the cubic axes a, b, c: 



a(b) 



= i(-a,f±V3a,r), 



(2) 



where and are Pauli matrices. The physical mean- 
ing of the T operators is to describe the fluctuations of 
exchange bonds due to orbital dynamics. On the cubic 
lattice, Eq. (0) can be rewritten in the equivalent form: 



J^-=4rfTf -2(Tr+r;) + l. 



(3) 



The main feature of this model is the strong inter- 
play between spin and orbital degrees of freedom which 
is suggested by the very form of Hamiltonian (||). In 
fact, the Kugel-Khomskii model contains rather nontriv- 
ial physics: That is, the classical Neel state in Eq. (H) 
{{SiSj) — —1/4) is infinitely degenerate due to the pres- 
ence of the orbital sector, and this extra degeneracy must 
be lifted by some mechanism. According to a recent pro- 
posal by Feiner, Oles, and ZaanenEl, this orbital frus- 
tration problem is likely solved by the formation of a 
RVB-type spin-singlet ground state. Technically, their 
suggestion is based on the observation that spin-orbit 
coupling results in a new, composite (simultaneous spin 
and orbital flip) excitation. This mode is found to be 
soft in certain directions in momentum space, leading to 
one-dimensional fluctuations, thus destroying the mag- 
netic order -completely. On the contrary, KhaliuUin and 
OudovenkoH have found the orbital flip excitation to be 
gapped, and concluded that the orbitally ordered quasi 
one-dimensional quantum Neel state is the proper low- 
temperature fixed point of model (||) . 

The purpose of our paper is twofold. First, we re- 
examine the solution of Feiner et alB and show that 
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FIG. 1. \3z^ — r^) orbital order which leads to weakly cou- 
pled AF spin chains (Jc — 4, J± = 1/4). As discussed in 
Ref. this type of orbital ordering provides the largest en- 
ergy gain due to quantum spin fluctuations. An orbital flip 
(indicated by an arrow) modulates the strength of the neigh- 
boring exchange bonds, breaking the c chain. In the classical 
Neel state, orbital excitations cost no energy. However, strong 
quasi one-dimensional spin fluctuations stabilize this struc- 
ture by opening the orbital gap through a order-from-disorder 
mechanism. 

the soft spin-orbital mode acquires in fact a firiite gap 
through Villain's order-from-disorder mechanisrr£3. Ba- 
sically, this gap is determined by the orbital excitation 
energy. Second, we calculate the spectral weight of the 
composite spin-orbital excitation in the dynamical spin 
structure factor, suggesting the possibility to measure the 
orbital excitation gap indirectly by a conventional neu- 
tron scattering experiment. 

To study model (§) we fqlJow the same scheme and 
method used by Feiner et alB. Namely, i) we start with 
the assumption of long-range staggered spin order and 
an uniform z-type orbital ordering (see Fig. |l|); ii) next 
we calculate the transverse spin-fluctuation spectrum by 
using the equation-of-motion method, and calculate the 
quantum corrections to the magnetic order parameter. 
We are motivated to use this method here for the fol- 
lowing reason: Employing a conventional diagrammatic 
technique it was argued in Ref. |^ that the orbital degener- 
acy of the classical state is removed by quantum spin fluc- 
tuations which generate orbital order and spontaneously 
break the cubic symmetry. The interesting question is 
then whether and how this order from disorder phenom- 
ena manifests itself in the equation of motion ansatz of 
Ref. |. 

The quantity to be calculated is the Green's function 
Gi.i' = {{S^\S^)), from which the dynamical spin struc- 
ture factor as well as the reduction of the spin order pa- 



rameter can be deduced. It obeys the following equation 
of motion: 

u;G.,,^{[St,S-]) + {{[S+,H]\Sr;)). (4) 

Within a conventional spin-wave approximation, this 
equation reads as 

(w - Am-j) Gi^i' = Sis' + 2mi ^ Gja' + ^m,; ^ Gj^i' 

jc j± 

+ (5) 

j± 

The following notations are introduced here: rui = 1 if 
site i belongs to sublattice A with up-spin orientation, 
and nii = — 1 otherwise. The summations over and 
j± are over nearest neighbors of site i in the c-chain 
and perpendicular directions, respectively. The factor 
\ = Jc + 2J^ = 9/2 and a„- = |i?f^.| - |i?f^-|. Dis- 
tinct from a Heisenberg model, a new Green's function 
Djii = {{Kjj\S^)) appears in the above equation. The 
operator KJj, — (JjS^, describes a simultaneous orbital- 
spin excitation which couples to a single spin flip. Fol- 
lowing Ref. H, we discarded in Eq. (||) the term Kj^ji, 
neglecting the correlation between orbital and spin exci- 
tations from different chains. 

The dynamics of the composite excitations K^^ — 
o'fS^ is described by the equation 

(- + 3... T i^j^ = -^(m. ± 1) E (xr, + f^sA 

i± 

+Lf, (6) 

where the last term is 

= E {^KpSj - ^KP^I ~ K T l)Kfi} ■ (7) 

la 

Here, site I is the nearest neighbor of site j, and 6Sj is 
the fluctuating part of the spin operator Sj. The 
term in Eqs. (||) and (0) accounts for the correlations 
within the c chain. If one neglects this term (as Feiner et 
aln did), one obtains from Eqs. (||) and ^ a soft mode 
mentioned above, which results finally in a breakdown 
of the assumed spin-ordered state. We do not accept 
this approximation. It is the intra-chain quantum fluc- 
tuations represented by the term that play a crucial 
role by making the orbital excitations (and conse- 
quently a^S'^) to acquire a flnite mass gap. The un- 
derlying physics can easily be observed in the following 
way: We write the equations of motion for the operators 
in Eq. (0) treating the bond (j7) exactly while consider- 
ing the remaining bonds of the c chain in a static Neel 
approximation. This results in 

- (m, ± 3)] KpS^ = (1 ± m,) ± KpS^) , 



{CO T mpSt = (1 T m.) (Xf, ± Kpsf) , 

[u - (3m. ± 1)] Kfi = 2KpSt - 2KpS] 

+(m,±l)X±. (8) 

Many terms in these equations in fact drop out due to 
I mil = 1, and Eqs. (||) together with Eq. (Q) lead to a 
very simple result. Namely, 

^tes-^Kl,, L-,, = -^KJ^, (9) 

and Lj^g = 0, L^^^^ = 0. Now, Eq. (^) can be written 
as follows (note K'' ^ K+ + R-): 

+^jK^r (10) 

The function = 8/(4mj + uj) is interpreted as a self- 
energy correction to the composite excitation cr^S'j' due 
to quantum fluctuations about the Neel state. We think 
that the large lu behavior of S is not quite reliable since 
short-time fluctuations of the environment of the bond 
(jl) were neglected in the above crude derivation. We 
therefore take the low-energy limit, Ej « 2mj — ^w, 
which is of main interest. Physically, the self-energy is 
mainly due to the orbital flip excitation (see Fig. P which 
acquires a finite energy as strong correlations within the 
c chains are explicitly taken into account. 

We now denote Gs — Gaa , Ds = Daa if the index s = 
1, and Gs — Gba, Ds = Dba if s = — 1. In momentum 
space, Eqs. (||) and ( |l0|) lead to 

(su; - X)Gs - XikG-s + VkD^s = 5s,i, 

[SLU-T- ^{SLU)] Ds - TJ^D^s + VkG-s = 0. (11) 

Here, E(lj) = 2 - lu/2, t = 3/2, jk = (8c, + 7±)/9, j± = 
(cx + Cy)/2, Vk = V^{cx - Cy)/2, and Cq = cosfc^. Eqs. 
( pi] ) can easily be solved to give two branches wi(fc) and 
0^2 (fc) for the spin excitations and their relative weights in 
a dynamic spin structure factor. The results are shown in 
Figs. H and p| for certain directions in momentum space. 

It should be noticed first that there is only one gapless 
mode at the F point. This is a plausible result since the 
original model has a continuous symmetry only in spin 
subspace; the second branch related to the breaking of 
the discrete orbital symmetry is expected to have a finite 
gap. We conclude that the soft mode discussed in Ref. 
^ is an artefact of the approximation which neglects the 
crucial effect of quantum spin fluctuations on the com- 
posite excitation a^S^. With regards to the low energy- 
momentum limit, the model (|^) behaves as a conventional 
Heisenberg system. However, the underlying anisotropic 
orbital ordering and the presence of orbital fluctuations 
strongly enhance quantum effects and reduce the order 
parameter. We have calculated (S^) = 0.23 which is con- 
sistent with Ref. ^. 




FIG. 2. Transverse spin excitation spectrum along the 
direction F - X(7r,0,0) - W(7r,f,0) - L(f , f , f ) - F - 
K(x, %^,0) in the Brillouin zone. To the right from L, the 
spin-orbit coupling vanishes and excitations are of either pure 
spin flip or of simultaneous spin-orbital flip character. Solid 
lines: present theory. Dashed line: spin-wave dispersion cal- 
culated with Jc = 4, J± = 1/4 and neglecting orbital fluctua- 
tions. Inset: on-site dynamic structure factor; the pseudogap 
seen about ci; ~ 2 is the manifestation of the orbital gap. 
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FIG. 3. Transverse spin excitation spectrum (lower panel) 
and the corresponding spectral weights in the dynamical 
structure factor (upper panel) along the F - X(7r, 0, 0) - 
N(7r, 0, tt) - F direction. Dashed line: spin-wave spectrum 
calculated neglecting spin-orbit coupling. 
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Second, at the magnetic zone boundaries the spin-orbit 
couphng results in quite visible deviations of the main 
magnon branch from that of an anisotropic Heisenberg 
model (see Figs. ||and^). This effect, which is due to the 
modulation of exchange bonds by orbital fluctuations, is 
expected to be a generic feature of spin-orbital models. 
In the present model, the zone boundary effect is related 
to the particular momentum dependence of the coupling 
constant 14 in Eq. ( pi] ) . Experimentally, anomalous zone 
boundary magnon softening has been observed in ferrp^ 
magnetic manganitesO. We have recently suggested^ 
that this effect originates from the modulation of double- 
exchange bonds by orbital fluctuations. 

Because of the finite hybridization of the excitations 
and S'^ (note that orbital pseudospin is not a con- 
served quantity), a conventional inelastic neutron scat- 
tering experiment might provide information on orbital 
excitation energies as well. In Fig. ^, we plot the relative 
spectral weight of two modes in the dynamical structure 
factor. A strong mixture of single spin S'^ and cr^5+ 
excitations and the "level repulsion" effect take place at 
energies where these excitations meet each other. In the 
present model, this results in an additional peak which 
has substantial weight at certain momenta. 

One comment is necessary at this point: The present 
approach (as well as the one of Rcf. ||) implicitly assumes 
the formation of a bound state of orbital and spin excita- 
tions (a tightly bound composite excitation). However, it 
might well be possible that the bound state decays into an 
independent orbital wave and magnons. In this case, one 
should expect instead of a well defined additional peak a 
softening and damping of magnons due to the coupling 
to the orbital- magnon excitation continuum. In fact, the 
latter picture is described in Ref. which treats the spin- 
orbit coupling perturbatively. In some low-dimensional 
spin-orbital models the bound state might existll3. The 
problem of a bound state versus particle-hole continuum 
in a Kugel-Khomskii model remains open at present. In 
either case, with regards to the low-energy behavior both 
approaches - the present one as well as that of Ref. |9| - 
are fully consistent with each other since the essence of 
the model, that is the orbital-order-out-of-spin-disorder 
mechanism, is equally taken into account. 

To conclude, we have studied a toy model providing a 
strong interplay between spin and orbital excitations. In 
the general context of transition metal oxides, this work 
suggests that the orbital fluctuations lead to a softening 
and damping of the zone-boundary magnons, and may 
even result in an additional structure in the spin-response 
fmiction of these compounds. 
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